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We have systematically investigated the magnetic properties and magnetocaloric effect (MCE) 
in i?Mn03 (-R=Dy, Tb, Ho and Yb) single crystals. Above a critical value of applied field {He), 
.RMnOs undergo a first-order antiferromagnetic (AFM) to ferromagnetic (FM) transition below 
the ordering temperature {T§) of ii"^^ moment and a second-order FM to paramagnetic (PM) 
transition above T^. Both H and T dependence of M shows that the system is highly anisotropic 
in the FM as well as PM states and, as a result, the magnetic entropy change (ASa/) is extremely 
sensitive to the direction of applied field and can be negative (normal MCE) or positive (inverse 
MCE). For hexagonal HoMnOa and YbMnOs systems, a very small inverse MCE is observed 
only for H parallel to c axis and it decreases with increasing H and crosses over to normal one 
above He- On the other hand, for orthorhombic DyMnOs and TbMnOs, though the inverse MCE 
disappears above along easy-axis of magnetization, it increases rapidly with H along hard-axis 
of magnetization for T<^T^. Except for YbMnOs, the values of ASm, relative cooling power and 
adiabatic temperature change along easy-axis of magnetization are quite large in the field-induced 
FM state for a moderate field strength. The large values of these parameters, together with 
negligible hysteresis, suggest that the multiferroic manganites could be potential materials for 
magnetic refrigeration in the low-temperature region. 



PACS numbers: 75.30.Sg,75.47.Lx,75.30.Kz,75.40.Cx 
I. INTRODUCTION 

In twenty-first century, energy efficient and environ- 
mentally friendly technology has received special atten- 
tion in order to combat the global warming phenomenon 
and energy crisis. Refrigeration based on the magne- 
tocaloric effect (MCE) has attracted much research in- 
terest because of its higher energy efficiency over the con- 
ventional vapor compression refrigeration and it does not 
use ozone-depleting chlorofluorocarbon as a refrigeranttli. 
Magnetocaloric effect describes the reversible change in 
temperature of a material under adiabatic condition pro- 
duced by the magnetic entropy ch ang e ASm due to the 
variation in applied magnetic fielcPSl. The main aim in 
this field is to search for new materials, which exhibit 
a large MCE and are capable of operating at different 
temperature ranges, depending on the intended applica- 
tions. Large MCE close to room temperature would be 
useful for domestic and several technological applications 
while large MCE in the low-temperature region is impor- 
tant for specific technological applications such as space 
science and liquefaction of hydrogen in fuel industrjUJ2l. 
The guidelines for the choice of an appropriate material 
are that it should have low heat capacity and exhibit a 
large entropy change at the ferromagnetic (FM) to para- 
magnetic (PM) transition or field-induced metamagnetic 
transition from antiferromagnetic (AFM) to FM states 
with a minimal hysteresis. 

In colossal magnetoresistive oxides, MCE has been 
extensively studied due to their wide variation in the 



Curie teinperature (Tc) and nature of FM-PM phase 
transitiorP"^. On the other hand, MCE in multifer- 
roic manganites i?Mn03 (i?=Tb to Yb) has not been 
investigated in detail^^. So far, most of the studies 
on these systems are concentrated on the magnetic and 
ferroelectric properties of Mn sublattice because their 
interplay has opened up a new dimension in basic re- 
search as well as technological application. The mag- 
netic structure determination is a prerequisite for under- 
standing the coupling between ferroelectricity and mag- 
netism in these materials. Several neutron and reso- 
nance x-ray studies reveal that RM11O3 undergo sequence 
of complicated magnetic phase transitions with decreas- 
ing temperatur^i21116]_ j^g^g q^^q been shown the rare- 
earth magnetic ordering plays a very important role in 
the magnetoelectric couplin^l^. According to the crys- 
tallographic structure, these multiferroic materials are 
divided into two classes. The compounds with larger 
rare-earth ion (i?=Tb,Dy) crystallize in orthorhombic 
structure whereas hexagonal structure is more stable for 
smaller ionic radius of R (i?=IIo to Lu,Y) . In orthorhom- 
bic compounds, where the magnetic frustration of Mn 
spin arises from competing exchange interactions, fer- 
roelectric and magnetic orderings appear at the same 
temperatur^i2]_ However, in hexagonal compounds, the 
magnetic frustration of Mn spin arises from the lattice ge- 
ometry since the triangular lattice is frustrated for AFM 
first-nearest-neighbor interaction and the ferroelectric or- 
der occurs at an elevated temperature (^900 K) which is 
well above the magnetic ordering temperature (~80 KJP^l. 
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Irrespective of the magnetic and ferroelectric properties 
of Mn sublattice, the AFM ordering of rare-earth mo- 
ments in multiferroic KMnO^, occurs at a relatively low 
temperatm'e (T^<8 K) and the magnetic structm'e is 
highly anisotropic with very weak interaction along one 
of the crystallographic axes. As a result, several mul- 
tiferroic manganites with high total angular momentum 
quantum number of rare-earth ion exhibit huge increase 
in magnetization at a moderate field strength (~2 T) 
which is close to the expected moment of R!^~^ ion. In 
this work, we present the field and temperature depen- 
dence of magnetic properties of rare-earth sublattice in 
i?Mn03 (i?=Dy, Tb, Ho and Yb) crystals. We observe 
that these materials (except R=Yh) exhibit giant MCE 
and large adiabatic temperature change and RCP due to 
the field-induced AFM-FM transition. This opens up a 
possibility for another viable technological application for 
multiferroic manganites namely, in magnetic cooling at 
low temperature. Furthermore, we have shown that one 
can clearly differentiate orthorhombic manganites from 
hexagonal ones because of their distinct magnetic and 
magnetocaloric properties related to rare-earth ion. 



II. EXPERIMENTAL TECHNIQUES AND 
SAMPLE PREPARATION 



Polycrystalline i?Mn03 (i?=Tb, Dy, Ho and Yb) sam- 
ples were prepared from stoichiometric mixture of -R2O3 
and Mn304 by solid-state reaction and the single crys- 
tals were grown from the polycrystalline rod by travelling 
solven t flo ating zone technique using an image furnace 
(NECj^Esl. Magnetic and heat capacity measurements 
were carried out employing superconducting quantum in- 
terference device magnetometer (Quantum Design) and 
physical properties measurement system (Quantum De- 
sign). The x-ray diffraction patterns of powdered sam- 
ple of single crystals reveal that these materials are sin- 
gle phase. The Rietveld profile refinement of diffrac- 
tion patterns shows that DyMn03 and TbMnOs have 
orthorhombic Pbnm structure whereas HoMn03 and 
YbMn03 exhibit hexagonal crystal structure with space 
group PGscm. Several magnetic parameters such as mag- 
netization (Af), susceptibility (x) and entropy (Sm) with 
field parallel to easy axis (-ff||e) are denoted as Me, Xe 
and Smc, respectively whereas the respective parameters 
with field parallel to hard axis {H\\h) are denoted as Mh, 
Xh and for both orthorhombic and hexagonal struc- 
tures. 
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FIG. 1. (a) Temperature dependence of dc susceptibil- 
ity along the easy-axis (6-axis) and hard-axis (a-axis) of 
DyMnOa. (b) Temperature dependence of xT^ along both 
axes. Solid lines are the linear fit to the data along the 
respective axis, (c) Thermal evolution of heat capacity of 
DyMnOa crystal at zero field. 



III. RESULTS AND DISCUSSIONS 

The isothermal magnetic entropy change IS.Sm can be 
calculated from H and T dependence of M using the 

Maxwell relation, ASm{T, H)^J {^) dH^ ^ . Since 



the magnetization measurements are performed at dis- 



crete intervals of T and H, ASm is numerically calcu- 
lated approximately using the expression, ASm{T, H) = 
Si ^-Ha, where M^+i and M^ are the experimen- 

tal values of magnetization measured with a field H^ at 
temperatures T^+i and T^, respectively. The characteris- 
tic parameter which determines the magnetic cooling effi- 
ciency of a magnetocaloric material is the relative cooling 
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power (RCP) and is defined as, 

RCP ^ I ASMdT 



200 



(1) 



where Ti and T2 are tire temperatures corresponding 
to both sides of the half-maximum value of ASi\i peak. 
The relative cooling power is the measure of the amount 
of heat transfer between the cold and hot reservoirs in 
an ideal refrigerator as a function of field. The adi- 
abatic temperature change (ATad), another important 
parameter related to MCE, can be calculated from the 
field-dependent magnetization and zero-field heat capac- 
ity data. The total entropy S{0,T) in absence of magnetic 
field is given by 



5(0, T) 



C(0,T) 



dT 



(2) 



and then S{H, T) may be evaluated by subtracting the 
corresponding ASm from 5(0, T). The isentropic tem- 
perature change between the entropy curves 5(0, T) and 
S{H,T) provides the value of ATaJyTf^. 



A. Orthorhombic DyMnOs and TbMnOa systems 

The temperature dependence of zero-field-cooled 
(ZFC) and field-cooled (FC) dc susceptibility x {^M/H) 
for H (=10 Oe) along the easy-axis (6 axis) and hard- 
axis (a axis) of magnetization are shown in Fig. 1(a). 
It is clear from the figure that DyMnOa is magnetically 
anisotropic and undergoes a PM to AFM transition be- 
low T^^^6.5 K due to the long-range ordering of the 
Qy3+ nioments. For further understanding the nature of 
magnetic interaction, we have plotted versus T [Fig. 
1(b)]. The linearity of x^^{T) over a wide temperature 
range suggests that susceptibility follows the Curie- Weiss 
(CW) law [x=G/{T — 6cw)]- From the high temperature 
linear fit, we have calculated the CW temperature Ocw^"^ 
K and the effective moment Pe//=11.8 /is for H\\e while 
the corresponding values are —66 K and 11.6 fiB for H\\h. 
The observed values of Peff are close to the theoreti- 
cally expected moment 11.7 /is, calculated using the two- 
sublattice model Pe// = y^(P^^p2 _^ [P^jjf. The large 

difference in O^.^ for two different crystallographic axes 
reflects the anisotropic nature of the exchange interac- 
tion. The small value of Ocw indicates that the magnetic 
interaction along the easy axis is very weak. On the other 
hand, the large and negative value of O^w indicate that 
the magnetic interaction along the hard axis is strong 
and antiferromagnetic in nature. For T^T^^ , the large 
difference in xC^) between two axes implies that the sys- 
tem remains highly anisotropic even in the PM state. Fig. 
1(c) shows the temperature dependence of specific heat 
in zero field. The anomalies in Cp at 38.2, 17.1 and 5.5 
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FIG. 2. Isothermal magnetization of DyMnOs as a function 
of magnetic field for different temperatures along the easy (a) 
and hard axes (b). Insets show the hysteresis at 2 K. Ratio 
of isothermal magnetization along the easy and hard axes 
{M^/Mh) as a function of H at some selective temperatures 
(c). 



K correspond to AFM transition of Mn moment {T^'") 
into a sinusoidal incommensurate phase, lock-in transi- 
tion (r||j") of Mn spin and AFM transition of Dy^+ ion, 
respectively. 

For elucidating the role of applied magnetic field on 
AFM ordering, we have measured H dependence of M 
in the vicinity of T^^ and beyond. The field is applied 
nearly parallel to the easy- and hard-axis of magneti- 
zation. For each isotherm, the magnetic field has been 
varied from to 8 T. Some representative plots of isother- 
mal field variation of M in the temperature range 2-42 
K are shown in Figs. 2(a) and (b), which depict a field- 
induced metamagnetic transition. From these plots, one 
can clearly differentiate the nature of H dependence of 



M below from that above Tt!". Below T^^, for both 
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FIG. 3. (a) Temperature dependence of magnetization along 
the easy-axis of DyMnOs for different magnetic fields. Inset 
shows the differential susceptibility (xde) versus temperature 
for different fields, (b) The Arrott plots with H along the 
easy axis, (c) The Arrott plots with H along the hard axis, 
(d) The H -T phase diagram of DyMnOs. 



LJ 1 I 1 I 1 I 1 I i_I 

10 20 30 40 
T(K) 

FIG. 4. Magnetic entropy change versus temperature for dif- 
ferent field change along (a) easy axis (b) hard axis. 



the axes, M increases slowly with H in the low-field re- 
gion followed by a sharp jump at a critical field He and 
then increases slowly with further increase of H. M does 
not show monotonic temperature dependence for H<Hc. 
This behavior is consistent with the field-induced transi- 
tion from AFM to FM state at H=Hc. We observe that 
He is slightly higher and the field-induced metamagnetic 
transition is less sharper for H\\h as compared to H\\e. 
Also, the value and nature of H dependence oi M ioi H\\h 
and H\\e are significantly different. This is more clearly 
refiected in the field dependence of the ratio Me/Mh, 
which shows a sharp and symmetric peak at around 1.5 
T [Fig. 2(c)]. As T increases, the peak becomes broad 
and shifts slowly towards higher H. Above 10 K, the 
peak changes into a shoulder-like feature. Another im- 
portant point to be mentioned here is that Mh{H) curve 
shows a weak feature at around 5 T below T^^; indicat- 
ing the presence of a second transition. Similar weak 
feature at high field has also been reported earlier in 
single crystals of DyMnOaP^. The insets of Figs. 2(a) 
and (b) display the five-segment M{H) loop at 2 K. We 
did not observe any hysteresis at low field. This is ex- 
pected because DyMnOa is AFM for H<1.5 T. However, 
polycrystalline samples with small grain size may show a 
weak hysteresis due to the surface ferromagnetisnl^^. As 
the size of crystallites is quite large in good quality sin- 
gle crystals, the contribution from surface magnetism is 
very small. The value of saturation magnetization (Mg) 
deduced at 2 K and 8 T is 8.8 hb for H\\e which is 88% 
of expected moment (10 /xb). In contrary to this, M 
does not show saturation-like behavior and its value is 



less than half of the expected moment for H\\h. The 
higher value of He, low magnetic moment and absence of 
saturation-like behavior are consistent with the large and 
negative value of 9cw for H\\h. The temperature depen- 
dence of M and differential susceptibility Xd {—dM/dH) 
are useful to understand the field-induced magnetic tran- 
sition. Figure 3(a) shows T dependence of M and Xd for 
H\\e. The peak position in Me{T) and Xde{T) curves 
corresponds to AFM transition temperature . T^^ 
decreases rapidly with increasing H and disappears at H^ 
in both the cases. Above H^, the nature of M{T) curve 
changes dramatically. M increases rapidly with decreas- 
ing T but saturates at low temperatures. The saturation 
region widens with the increase of field strength. These 
behavior suggest that DyMnOa undergoes field-induced 
AFM-FM transition below T^^ and PM-FM transition 
above . Normally, the field- induced order-order tran- 
sition is first order in nature. For understanding the 
nature of AFM-FM and FM-PM phase transitions, we 
have transformed the Mf,{H) data into Arrott plots as 
shown in Fig. 3(b;P. The slope of H/M versus JVP 
curve is useful to determine the order of both temper- 
ature and field driven magnetic phase transition. The 
negative slope of the Arrott plot often indicates a first- 
order nature of the transition, while the positive slope 
implies a second-order transitiorP^l. In the present case, 
the negative slope below T^^ for H<1.8 T indicates that 
the field-induced AFM-FM transition is first-order in na- 
ture while the positive slope of the high-field data and 
their linear extrapolation to H=0 at a non-zero positive 
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value of M for T^^<T<2Q K suggest the second-order 
nature of the PM-FM transition. The Arrott plots are 
also done for Mh{H) data [Fig. 3(c)]. Though the Ar- 
rott plots for H along two different directions are quali- 
tatively similar below T^^, the nature of high- field data 
above T^^ in two cases are very different. The first-order 
nature of AFM-FM transition is also evident from the 
low-field Mh data. However, the linear extrapolation of 
high-field Mh data to H—Q does not reveal any non-zero 
and positive value of M for T>T^^ . This suggests that 
the PM-FM phase boundary at T^^ is extremely sharp 
when the field is applied along the hard-axis of mag- 
netization. The temperature dependence of He, deter- 
mined from the maximum in dM{H)/dH (for T<T^^) 
and minimum in dM(r)/dr (for T>T^y) curves and the 
Arrott plots, is summarized in the (H—T) phase diagram 
[Fig. 3(d)]. In order to test whether DyMnOa is a suit- 
able candidate for magnetic refrigeration, we have calcu- 
lated the isothermal magnetic entropy change from the 
M{H) curves [Figs. 2(a) and (b)] using the Maxwell rela- 
tion. Figure 4 presents the thermal distribution of A5j\/ 
for field variation up to 8 T for H parallel to easy and 
hard axes. In both the cases, AS" a/ is negative above T^^ 
and the magnitude of -ASm{T) at maximum (A5'™^) 
increases with field. AS'JJj;'*^ is as high as 5.52 J mol~^ 
K"-'^ at 8 T for H\\e. However, the H dependence of 
ASm below T^^ are quite different in two cases. ASme 
is positive (inverse MCE) below T^^ only for small field 
changes AH<Hc (=1.5 T). On the other hand, ASmh 
remains positive up to 8 T for r<3 K. Moreover, the 
steep decrease in ASMhiT) on the low-temperature side 
of the maximum suggests that DyMnOa may exhibit a 
large inverse MCE below 2 K. It is worth noting that 
the position of the maximum in -ASm{T) curve shifts 
slowly towards higher temperature with increasing field 
for H\\e while it is insensitive to field for H\\h. This be- 
havior is consistent with the observed {H — T) phase di- 
agram. For H\\e, the positive value of ASmc below T^^ 
initially increases and then decreases with the increase 
of H and eventually becomes negative above 1.5 T. Such 
type of H dependence of ASm is quite common in anti- 
ferromagnetic systems where the field-induced AFM-FM 
transition occurs^. The initial increase in ASm with H 
below He is due to the field-induced magnetic disorder- 
ing. When an external magnetic field is applied along the 
easy axis, the magnetic moment fluctuation is enhanced 
in one of the two AFM sublattices which is antiparallel 
to H. With the increase of H , more and more spins in 
the antiparallel sublattice orient along the field direction. 
This, in turn, increases the spin disordering and it con- 
tinues up to a certain field below He- As the system be- 
comes ferromagnetic for H>Hc, the majority of spins in 
the antiparallel sublattice orient along the field direction 
and, as a consequence, ASm becomes negative. The rel- 
ative cooling power is evaluated to determine the cooling 
efficiency of DyMnOs crystal as a magnetocaloric mate- 
rial. RCP is quite large (155 J mol~^) for a field change 
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FIG. 5. (a) Temperature dependence of FC and ZFC dc sus- 
ceptibility curves for HoMnOa along easy axis (a axis) and 
hard axis (c axis), (b) Inverse susceptibility versus tempera- 
ture measured along easy and hard axis. Solid line indicates 
the Curie-Weiss fit in both directions, (c) Temperature de- 
pendence of specific heat for HoMnOa crystal at zero magnetic 
field. 



of 8 T for H\\e. We have also estimated the adiabatic 
temperature difference from the isentropic curves using 
the zero-field specific heat data in Eq. (2). ATad is found 
to be as high as 11.5 K for a field change of 8 T for -ffJIe. 
This value of ATad is appreci ably larger than that ob- 
served in perovskite manganite^^^. Similar to DyMnOa, 
the temperature and field dependence of magnetization 
with H along easy-axis (a axis) and hard-axis (b axis) of 
TbMnOa crystal has been measured. We observe that 
the T and H dependence of magnetization for TbMnOa 
crystal is similar to earlier reports^ Unlike DyMnOa, 
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FIG. 6. Figure 6. Field dependence of isothermal magnetiza- 
tion for HoMnOa with H parallel to hard axis (a) and easy 
axis (b). Me/Mh versus H plot at some selective tempera- 
tures (c). The Arrott plots of the magnetization isotherms 
along hard axis (d) and easy axis (e). 



the inverse MCE in TbMnOa is appreciably large and 
increases with H along h axis. This suggests that the 
AFM ground state in TbMnOs is quite stable against 
applied field along the hard axis of magnetization. In 
spite of this difference, the nature of field-induced AFM- 
FM and PM-FM transitions and [H — T) phase diagrams 
for H\\e as well as for H\\h are qualitatively similar in 
two systems. Also, we observe that the values of magne- 
tocaloric parameters A5]J}''^, RCP and AT^d in DyMnOg 
and TbMnOs are comparable in magnitude. Therefore, 
both DyMnOs and TbMnOa satisfy the major important 
criteria for magnetic refrigeration at low temperatures. 



B. Hexagonal HoMnOs and YbMnOa systems 

The temperature dependence of low-field {H—A Oe) 
susceptibility along the two principal axes of hexagonal 
unit cell of HoMnOa is shown in Fig. 5(a). The anomaly 
at 4.5 K in ZFC cycles is due to the AFM ordering of Ho 
moments whereas the weak bifurcation between FC and 
ZFC curves below ^^40 K corresponds to the reorientation 
of Mn moment in the basal plane perpendicular to the ini- 
tial direction. No anomaly is observed at the AFM order- 
ing temperature of the Mn moments. Figure shows that 
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FIG. 7. Figure 7. Temperature variation of magnetic entropy 
change for different field change along (a) easy axis (b) hard 
axis for HoMnOs. 



the values of x ^-t any temperature are significantly differ- 
ent along two crystallographic directions; reflecting the 
anisotropic nature of the magnetic structure of HoMnOa. 
For the quantitative estimation of anisotropy in magnetic 
parameters, we have presented x~^iT) for H along easy 
axis (a axis) and hard axis (c axis) [Fig. 5(b)]. Both the 
curves show that x foUows the CW behavior over a wide 
range of T. From the linear part of X~'^{T) for H\\h, the 
fitted values of P^f / and Ocw are found to be 10.24 hb and 
— 116 K, respectively while the corresponding values are 
11.2 /is and — 5 K for iJ||e which are consistent with the 
previously reported result^^HEfi] Significant anisotropy in 
magnetic interaction is evident from the large difference 
in the values of Peff and Ocw for two principal axes of 
the unit cell. Unlike susceptibility, the temperature vari- 
ation of zero-field heat capacity is showing three distinct 
transitions [Fig. 5(c)]. The peak at 4.6 K indicates the 
AFM ordering of Ho moment as evidenced from magnetic 
measurement. The small peak at ^39 K is due to the 
Mn moment reorientation and the third A-type anomaly 
originates from the AFM ordering of the Mn'^+ magnetic 
moments. Figures 6(a) and (b) present the field depen- 
dence of isothermal magnetization for HoMnOa with field 
parallel to c and a axes, respectively. We would like to 
mention that M (at 3 K and 5 T) is 20-25% larger for 
the present sample than the reported value on HoMnOa 
single crystal grown by the flux methotpl. The M{H) 
isotherms for H\\h indicate a field-induced metamagnetic 
transition. At low fields, the dependence of M on is ap- 
proximately linear and the slope oi M{H) curve decreases 
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FIG. 8. Figure 8. Isothermal magnetization of YbMnOa as a 
function of magnetic field for different temperatures with H 
along the a axis (a) and c axis (b). 



above a critical field but no saturation is attained up to 
8 T. To investigate the reversibility of the field-induced 
magnetic transition, we have measured the five-segment 
M{H) loop at 2 K and observed no significant hystere- 
sis [Inset of Fig. 6(a)]. The nature of M{H) curves for 
H\\e differs from those for H\\h. For H\\e, M{H) curve 
displays no abrupt change but increases smoothly with 
H and the value of M is significantly larger as compared 
to that for H\\h. The anisotropy in the field-induced 
FM state is apparent from the Mf,/Mh versus H plot 
for different T [Fig. 6(c)]. The field and temperature 
dependence of A/g/M/i for HoMnOs may be compared 
with that for orthorhombic DyMnOs. One can clearly 
see that there are several important differences between 
the two systems. For example, in DyMnOa, the peak in 
Me/Mh versus H plot is very sharp at low temperature. 
On the other hand, the maximum for HoMnOs is quite 
broad and its height is almost half of that for DyMnOs. 
Also, the nature of T dependence of the peak are quite 
different in two cases. With increasing T, though the 
peak broadens, the value of Me/Mh at peak increases 
slowly in HoMnOa while the peak decreases and broad- 
ens rapidly in DyMnOa. Figures 6(d) and (e) show the 
Arrott plots for field along both the axes. The Arrott 
plots show that the field-induced AFM-FM transition is 
first order while PM-FM transition is second order and 
the PM-FM phase boundary is quite sharp for field along 
hard-axis of magnetization [Fig. 6(d)]. However, closer 
inspection shows that the nature of the Arrott plots for 
H\\e [Fig. 6(e)] differs significantly from that for H\\h 
and also from that for orthorhombic systems. For H\\e, 
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FIG. 9. Figure 9. The Arrott plots of the magnetization 
isotherms for YbMnOa with field direction parallel to a axis 
(a) and c axis (b). 



H/M versus Af^ curves do not show negative slope in 
low-field region as expected for the first-order transition. 
This is due to the smooth increase of M with H [Fig. 
6(b)]. For H\\h, the H/M versus curves for HoMnOa 
are almost linear in the high-field region both above and 
below T^°. Normally, the H/M versus M curves are 
linear and parallel for mean-field like FM-PM transition 
otherwise a upward curvature is exhibited as in the case 
of H\\e. A small deviation from parallelism may occur 
due to the weak field dependence of Tc in the high-field 
region. Thus it appears that the system may belong to 
different universality class depending on the direction of 
applied magnetic field. For HoMnOa, the isothermal en- 
tropy change for different field variations was estimated. 
Figures 7(a) and (b) depict the nature of ASm{T) curves 
ioi H\\e and H\\h, respectively. In both the cases, a broad 
maximum is observed. Apart from the values, there are 
several important differences in the nature of H depen- 
dence of ASm for two crystallographic axes. With the in- 
crease of magnetic field, the maximum in ASMiT) shifts 
to higher temperature for H\\e while the position of the 
maximum remains pinned at ~5 K for H\\h. For H\\e, 
ASm is negative down to 2 K even at very low fields, 
i.e., inverse MCE is absent in the AFM state. However, 
for H\\h, ASm is positive well below at low fields 
{H<Hc). Unhke DyMnOs, inverse MCE in HoMnOa is 
very small. One can see that AS^"^ is large for a moder- 
ate field change along easy axis. The values of other mag- 
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FIG. 10. The magnetic entropy change in YbMnOs for a field 
change from 0-1 to 0-8 T along the a axis (a) and c axis (b). 



netocaloric parameters such as RCP and ATad are also 
large for HoMnOa single crystals. We observe RCP=144 
J mol-i and Ar<j<j=12.5 K for the field change of 8 T 
along easy axis. 

We now briefly discuss the magnetic and magne- 
tocaloric properties of YbMnOa crystal. Similar to other 
multiferroic manganites, magnetic properties of rare- 
earth sublattice in YbMnOs is sensitive to the direction 
of applied field with respect to crystallographic axis. Fig- 
ure 8 shows the field dependence of magnetization for 
YbMnOa at different temperatures. One can see that 
the nature of magnetic response of the crystal with field 
along a and c axes are very different. Unlike HoMnOs, 
M is slightly larger along the c axis as compared to basal 
plane and two field-induced transitions are clearly visi- 
ble. However, M along a axis increases smoothly with 
H similar to that in HoMnOa crystal. The transition 
at high-field (^-^3 T) is extremely sharp (step-like) below 
(=3.6 K), shifts progressively towards higher field 
with increasing T and clearly visible up to 14 K [Fig. 
8(b)]. It is also interesting to note that M does not 
show saturation-like behavior up to 8 T and the value 
M at 2 K and 8 T for both the axes are comparable and 
about 37-39% of the expected moment. This behavior is 
quite different from other multiferroic manganites where 
M along two axes are significantly different and the ob- 
served high-field value of M for H\\e is close to expected 
moment. This suggests that YbMnOs is magnetically 
less anisotropic and the AFM interaction along both the 
axes is strong. Indeed from the magnetization measure- 
ments in the paramagnetic state, we observe that 6^^ is 
large and negative along both the crystallographic axes 



and their values are 195 and 225 K within the basal plane 
and along c axis, respectively. In order to determine the 
order of field-induced magnetic transition, we have also 
done the Arrott plots of magnetization data of YbMnOa 
[Fig. 9]. The nature of the Arrott plots for H along a 
axis is qualitatively similar to those for HoMnOa crystal 
but for H along c axis the low-field peak is quite promi- 
nent and the change in slope from positive to negative 
is abrupt, marked by a much sharper triangular-shaped 
minimum at He- However, unlike HoMnOs, the linear 
extrapolation of high-field magnetization data to H=Q 
reveals a non-zero positive value of M up to 2Tjj^ for H 
along c axis. Figure 10 illustrates the temperature depen- 
dence of A5'j\/ in a magnetic field change up to 8 T. A5m 
reaches a negative maximum value ~2.3 J mol~^ K"-'^ for 
AH=8 T. Though the magnitude of AS'm is smaller for 
YbMnOa due to the smaller total angular momentum 
quantum number, the overall nature of T and H depen- 




FIG. 11. The maximum of magnetic entropy change (a), max- 
imum adiabatic temperature change (b) and relative cooling 
power (c) as a function of field for i^MnOs crystals with H 
parallel to easy axis. 
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dcncc of ASm is similar to HoMnOa. We observe that 
the value of AT^d is significantly large (^15 K) in spite of 
its smaller ASm and RCP (26 J mol^^). In order to com- 
pare and contrast the nature of MCE among multiferroic 
manganites, the variation of AS^""=, AT™"=^ and RCP for 
H\\e have been plotted in Fig. 11 which shows that these 
parameters increase monotonically with AH. Though the 
values of AS^"^, ATJ^j*^ and RCP depend on the system, 
the nature of their field dependence is qualitatively sim- 
ilar. It is clear from the figures that among the four sys- 
tems we have studied, DyMnOs, TbMnOs and HoMnOa 
are most suitable as low temperature refrigerant due to 
the high values of magnetocaloric parameters at a mod- 
erate field change. 

IV. CONCLUSIONS 

In conclusion, the detailed analysis of magnetization 
data shows that the magnetic interaction within the 
rare-earth sublattice in multiferroics -RMnOs is highly 
anisotropic (except i?=Yb). Above a critical field He, 
i?Mn03 undergo a field-induced first-order metamagnctic 
transition from AFM to FM state and a second-order 
PM-FM transition along with huge magnetic entropy 
change. Depending on T and the direction of applied 
field, ASm can be negative or positive, i.e., i?Mn03 ex- 
hibit both normal and inverse MCE. Except at low tem- 



peratures well below Tjj, ASm is negative and large in 
the field-induced FM state. For orthorhombic DyMnOa 
and TbMnOa, the inverse MCE is small, decreases with 
H and eventually crosses over to normal one above He for 
H\\e while it is large and increases with H for H\\h. In 
contrary to this, hexagonal HoMnOa and YbMnOs may 
show small inverse MCE only below H^. and at very low 
temperature for H\\c. Both orthorhombic and hexagonal 
multiferroics, except YbMnO.j. exhibit giant MCE, and 
large adiabatic temperature change and relative cooling 
power for H\\e. The large values of these parameters, 
negligible hysteresis and highly insulating nature suggest 
that the multiferroic manganites could be potential ma- 
terials for magnetic refrigeration in the low-tcmporature 
region. The present results also show that orthorhombic 
and hexagonal multiferroic manganites can be differenti- 
ated based on the magnetic and magnetocaloric proper- 
ties of rare-earth sublattice. 
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